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A full treatment for the scattering of an arbitrary number of bosons through a Bell multiport beam 
splitter is presented that includes all possible output arrangements. Due to exchange symmetry, the 
event statistics differs dramatically from the classical case in which the realization probabilities are 
given by combinatorics. A law for the suppression of output configurations is derived and shown to 
apply for the majority of all possible arrangements. Such multiparticle interference effects dominate 
at the level of single transition amplitudes, while a generic bosonic signature can be observed when 
the average number of occupied ports or the typical number of particles per port is considered. The 
results allow to classify in a common approach several recent experiments and theoretical studies 
and disclose many accessible quantum statistical effects involving many particles. 

I. INTRODUCTION 

Non-interacting, distinguishable particles exhibit independent and therefore uncorrelated behavior. Due to the 
bosonic or fermionic nature of identical particles, however, such statement is no longer true for indistinguishable par- 
ticles, even if no interaction takes place. For example, the bosonic nature of photons is impressively demonstrated by 
their statistical behavior in a Hong-Ou-Mandel (HOM) setup [1]. Here, two identical photons are sent simultaneously 
(within their coherence time) through the two input ports of a balanced beam splitter. Due to the lack of interaction 
between the photons, one would not expect any correlations in the number of photons measured at both output ports. 
For fully indistinguishable photons, however, the particles always leave the setup together, and never exit at different 
ports. 

Such synchronization of two non-interacting particles has lead to many applications in quantum information sciences. 
The visibility of the HOM-dip quantifies the indistinguishability of two photons [2 . Thereby, the quality of single- 
photon sources can be tested [3]. The maximally entangled |^~) Bell-state (in any degree of freedom carried by the 
photons) can be detected or created, since it leads to an unambiguous signature in the setup [1]. This projection onto 
an entangled state can be applied in entanglement swapping protocols [5] and quantum metrology [5]. 

It is therefore of great interest to generalize the HOM setup for more than two photons and more than two input 
or output ports, i.e. to n particles that are scattered in a setup with n input and output ports. This would allow 
applications such as entanglement swapping or entanglement detection for many particles and the experimentally 
controlled transition from indistinguishability to distinguishability for many identical particles [71 15] . 

While a comprehensive understanding of this scattering scenario is not yet available due to the complexity of 
the problem and the prohibitive scaling of the number of output states, several steps have been undertaken in this 
direction. The measurement of the enhancement of events with all particles in one port - bunching events, was realized 
experimentally [HI HH] j a prediction for the suppression of coincident events for a specially designed biased setup with 
three particles and three input ports was presented [TT] . The case of a Bell multiport beam splitter [T2 US] which 
redistributes n incoming particles to n ports in an unbiased way was discussed in [14 , where it was shown that 
coincident events are suppressed when n is even. 

In this contribution, we extend our recent results [15j on the characterization of the probabilities of all possible 
output events of the Bell multiport beam splitter when n particles are prepared in the n input ports. Such treatment 
enables a general understanding of multiparticle interference effects, as well as on the average behavior of bosons. It 
hence unifies previous experimental and theoretical work on multiport beam splitters, and opens up new perspectives 
for the experimental verification and exploitation of bosonic multiparticle behavior. 
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II. FORMALISM 



A. Setup and notation 



We consider a scattering scenario in which n particles are initially evenly distributed among n input modes. They 
are scattered by the multiport beam splitter and exit among n output ports. The probability for each particle to exit 
at any port is the same, i.e. 1/n. Such setup can be realized for photons via a simple beam splitter in the two-photon 
HOM case, as illustrated in Fig. ([T^). A pyramidal combination of beam splitters with different reflection/transmission 
rates yields the generalization for n ports and particles, such setup is shown in Fig. ([ijj), for n — 5. 

(a) (b) 




FIG. 1: (a) Two-photon HOM scattering setup with two photons impinging on a balanced beam splitter, (b) GeneraUzation: 
multiport beam splitter, here with five incoming and outgoing modes. 

We denote arrangements of n particles in the n modes by a vector s = (si, S2, . . . s„), with Sk the number of particles 
in the output mode fc, and X]r=i ~ ^- Alternatively, we define the port assignment vector d of length n with entries 
that specify each particle's output port. It is constructed by concatenating sj times the port number j: 

e.g., for the arrangement s — (2,1,0,2,0), we find d{s) = (1,1,2,4,4). 



B. Distinguishable particles 

For distinguishable particles, no many-particle interference takes place, and the probability for a certain arrangement 
s is given by simple combinatorics: 

1 n\ 
classic) = n^n;=iSj! 

Due to the lack of interference phenomena, we call this situation "classical" . In accordance to our intuition, probabili- 
ties are summed instead of amplitudes. Two classes of events are especially interesting due to their extremal character. 
Coincident events, i.e. Sc — (1, 1, . . . 1), are realized with probability n!/n". Bunching events, with all particles at 
one output mode fc, correspond to Sfe = n and thus to Sb = (0, 0, .., n, ..0). They are realized with probability 1/n" 
and suppressed by a factor of n\ with respect to the coincident events. For large n, both events are highly unlikely, 
extreme cases. 



C. Indistinguishable bosons 



We reformulate the scattering scenario for identical particles, in second quantization. Applications of our study are 
feasible with today's optical technologies [H], therefore, we focus on bosons. The initial state with one particle in 
each input port reads j^*) = Y[7=i '^I I'-') • The single-particle unitary evolution induced by the scattering setup acts 
on all particles independently and maps the input port creation operators a| to output creation operators 6| via a 
unitary matrix U [16 , such that 6^ = X]fc=i ^jkO^k- The unbiased Bell multiport beam splitter under consideration 
here corresponds to the unitary operation given by the Fourier matrix, defined for any dimension n by 

P^O-i)(fc-i) 

U,k^- j= . (3) 
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The possible states with fixed particle number per port after the scattering process read 



(4) 



The transition probability to a specific output arrangement s can be written with the help of the port assignment 
vector (Eq. [T]) as 



(5) 



where P„ denotes the set of all permutations of {1, .., n}. This coherent sum over n\ terms expresses the interference 
that occurs between all many-particle amplitudes that lead to the same output state. 



D. Equivalence classes 

In order to discuss the behavior of the scattering system, it is necessary to identify classes of final states that occur 
with equal probability, within the quantum and the classical case. In the latter case, the realization probability of any 
arrangement s, Eq. ([2|, remains invariant under permutation of the output ports Sk- Hence we can define classical 
equivalence classes which identify arrangements related to each other by permutation. Ultimately, however, all final 
events s have to be considered as inequivalent. In the case that the scattering matrix is a Fourier matrix such as given 
m Eq. ^, some symmetry properties allow to reduce considerably the number of equivalence classes. One indeed 
finds that the amplitude ([s]) is invariant under cyclic and anticyclic permutations. This allows us to define a quantum 
equivalence relation between arrangements, and N^y^g^^^ associated quantum equivalence classes. 

The number of classical equivalence classes corresponds to the partition number, i.e. the number of possibilities 
to write an integer as sum of positive integers, while the total number of inequivalent events grows much faster, it is 
given by N^q^^i = 2{n\)^ ■ ^'-'^ comparison, the number of equivalence classes are given in Table I 



III. EVENT-SUPPRESSION LAW 



In general, the evaluation of the transition probabilities in Eq. ([5| is a difficult task and cannot be performed in 
polynomial time with n [17) . It is, however, possible to exploit the symmetry of the matrix U to formulate a powerful 
law which predicts the suppression of final events. Indeed, since only n-th roots of unity appear in the Fourier matrix 
(Eq. [3]), also every term of the n! summands in Eq. ([s]) can be written as such. Thereby, Eq. ([s]) turns into 

n-l 

(vI/|<i>(sO) = ^Cfce'^^ (6) 

k=0 

where the are natural numbers which give the cardinality of the following sets, defined in analogy to |18] . 

©ri,s(CT) = ^ d{s)ia{l) = r mod n \ , (7) 
1=1 J 

with Cr = |ur(s)|. The sum corresponds to the position of the barycenter of the set of points {cfce'^'^|fc G {1, 
in the complex plane. We set Q — mod ^li^i ^'^d define an operation 7 which acts on permutations such 

that 7((T)(fc) = (T(fc) + 1 mod n. We find that On, (7(0-)) = 6n,s (cr) + Q- Thus, if Q 7^ 0, the repeated appHcation of 
7 gives us a bijection between all pairs of Ur+a-Q, for a G {0, 1, .., n — 1}. Hence, we find 

Vr e {0,..,n- l},Va G N : Cr+a-Q = c^- (8) 

If Q 7^ 0, the set of points {cfce*^'^|fc S {1, .., n}} describes several interlaced polygons centered at the origin, ensuring 
that the sum vanishes, hence the process with the final state s is suppressed in this case. Thus, without knowing 
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the values of the individual Ck, and only by symmetry properties, it is possible to predict that the total sum (Eq. [5]) 
vanishes. This observation allows us to formulate: 

n 

Q(s) :=Mod(^d,(sO,?^) ^0 ^ (*|$(.^)=0. (9) 
1=1 

The law can be applied on any final state in an efficient way; consider, e.g., n = 6 and si = (2,1,2,1,0,0). 
The port assignment vector reads (1, 1, 2, 3, 3, 4), and one finds Q{si) = 2, and this event is hence strictly suppressed. 
Unexpectedly though, the event 5*2 = (0, 1, 2, 0, 2, 1), which is obtained from si by simple permutation, gives Q(s2) = 
due to the different port assignment vector (2,3,3,5,5,6). It is actually enhanced by a factor larger than seven as 
compared to the classical event probability (also see Table 



A. Suppressed arrangements 

It is possible to estimate the number of suppressed arrangements predicted with the help of (|9| by a simple argument. 
Since the number of arrangement -Equant much larger than n, we can assume that the Q{s) are uniformly distributed 
in the interval [0, . . . ,n — 1] for the ensemble of events s. Then the probability to find a suppressed arrangement is 
given by the weight of nonvanishing values of Q{s}, i.e., by [n — l)/n ~ 1 — 1/n. This estimate is also numerically 
confirmed in the values shown in Table H] 



B. Application of the suppression law 

For n — 2. .6, we list the unsuppressed arrangements in Table [IT] together with their quantum enhancement, i.e., 
the ratio of quantum-to-classical event probability. From the results, it is apparent that the behavior of the system 
becomes more extreme, the more particles are involved: the enhancement factor is bounded from above by rt!, a 
value that is reached for the enhancement for bunching events with s — (n, 0,..,0). For the 50 quantum equivalent 
arrangements that exist for n ~ 6, we show the classical and quantum probabilities in Fig. [2j Two arrangements are 
suppressed although they do not fulfill the requirements of the law: s = (0, 1, 2, 1, 0, 2) and s = (0, 1, 1, 2, 1, 1). In Fig. 
[3j we show the values of the corresponding in the complex plane. One can easily see that while the values do not 
lie on polygons, the sum of all contributions still vanishes. Such situations are exceptional, as can be seen from Table 

IB 
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FIG. 2: Quantum (blue circles) and classical (red squares) probability for the realization of the 50 different quantum equivalence 
classes for n = 6 particles, sorted by their classical realization probability. Note that most states are fully suppressed in the 
quantum case, only a few are highly enhanced with respect to the classical case. 



IV. BOSONIC BEHAVIOR 



As we have seen in the last section, the implications of Eq. ([9]) on the realization probabilities of single events are 
important for the overall behavior of the system: most events are totally suppressed, while only few remain which 
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TABLE I: Total number of events (-^total)' classical equivalence classes 
(-'^class)' quantum equivalence classes (A'^quantum)' classes that therewithin 
fulfill the law PI (A'jg^^), and suppressed classes which are not predicted by 
Eq. [9] (iVsuppF 
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TABLE IL Nonsuppressed output states, together 
with the corresponding quantum enhancement, i.e., 
the ratio of quantum to classical event probability. 



are highly enhanced. Intuitively, one would expect that events with many particles in one port are generally favored 
by bosons. Indeed, bunching events are always enhanced by a factor of n\ with respect to the classical case. The 
number of particles in one port or the number of occupied ports does, however, turn out not to be a good indicator 
for the enhancement or the suppression of a certain event. For example, events of the type s = (n — 1, 1, 0, .., 0) could 
be expected to be enhanced due to the bosonic nature of the particles, while they actually turn out to be strictly 
suppressed, for all n. Thus, at the level of the event probabilities of single arrangements, interference effects dominate, 
and the bosonic nature of the particles is not apparent at all. 

Such general bosonic behavior is recovered when a coarse-grained grouping of many final arrangements in larger 
classes is performed. Such classes can be characterized, e.g., by the number of occupied ports k, by the number m 
of particles in one port, or by the classical equivalence classes. The event probability for such a class is given by the 
sum of the probabilities of the single events that pertain to the class. 

When performing such average, we expect that interference effects disappear while the bosonic enhancement of states 
with many particles in one port persists. This can be also seen in our formalism: according to ([5|, the probabilities 
Pqm{s) are given in terms of a sum over permutations of scattering amplitudes^ i.e., over complex numbers of equal 
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FIG. 3: Illustration of the corresponding Ck in the complex plane for the process with s — (0, 1, 2, 1, 0, 2) (left hand side) and 
s — (0, 1, 1, 2, 1, 1) (right hand side). The sum of the complex numbers vanishes while the points do not lie all on circles. 



modulus (products of matrix-elements of Ujk)- Since these numbers typically have different phases, they tend to 
add up destructively. However, all Sjl permutations a that interchange the Sj particles that exit in port j leave the 
scattering amplitudes invariant, so that Sj\ terms in the sum have equal phases and add up constructively. This 
motivates the following approximation for the transition probability ([s]): 



approx 





^class(^ 






^class(^ 



(10) 



We show the probability distribution for the number of occupied ports, for the classical calculation (Eq. [2]), for the 
bosonic quantum case (Eq. [5|, and for our approximation (Eq. 10), for n — 14, in Fig. |4] As expected, bosons always 
tend to occupy less output ports than in the classical case. This behavior is persistent for any n. Furthermore, Fig. [4] 
shows that the approximation (Eq. |10p predicts the actual outcome very well for most fc, and only fails for events with 
almost all or almost no sites occupied. This is easily understood, since, for very small or very large fc, few distinct 
equivalence classes contribute to these event groups. Then, again interference dominates the event probability, rather 
than bosonic behavior. 
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FIG. 4: Event probability for a given number of occupied ports, for n = 14. Red squares denote classical combinatorics, blue 
triangles the quantum calculation, and black circles our (bosonic) estimate for the quantum result. The inset shows the same 
distribution on a linear scale. Events with 14 occupied ports are strictly suppressed in the quantum case. 



Also the event probability for a given number of particles in one single port is well described by our estimate (Eq. 
10). For 14 particles, the probability distribution is shown in Fig. [sj Again, we see a dramatic difference between the 



classical and quantum case, especially for the probability to find a large number of particles in one port. 

A further application of the bosonic approximation can be performed when we group the 718146 events according 
to their 135 classical equivalence classes. The resulting probabilities are shown in Fig. [6j While this grouping is much 
more fine than in the previous two examples, the difference between the classical and quantum case is still very well 
pronounced and reproduced by the estimate. 
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FIG. 5: Probability to find exactly k particles (liorizontal axis) in one port, for n = 14. Red squares denote the classical, blue 
triangles the quantum calculation, and black circles the estimate ( |10[ ). The inset shows the distribution on a linear scale for 
small k. Events with 13 particles in one port are totally suppressed for bosons. 
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FIG. 6: Quantum (Blue triangles) and classical probability (Red rectangles) and estimate (Black circles) for the realization of 
events grouped according to the 135 classical equivalence classes. The classical equivalence classes are ordered according the 
classical realization probability given by Eq. (|2|. The classical probability is therefore - by construction - a monotonically 
increasing function. The inset shows a detail of the probability distribution in linear scale. 



V. CONCLUSIONS AND OUTLOOK 



Synchronization of non-interacting particles might seem contradictory by construction. However, it turns out to 
be possible due to the exploitation of quantum statistical effects stemming from the indistinguishability of particles. 
We generalized the most prominent example for such a behavior, the HOM effect, to n particles and n ports on two 
different levels: interference effects inhibit the realization of most possible events for single transition amplitudes, 
while general statistical characteristics with smooth bosonic behavior emerge which are efficiently approximated by 
Eq. ( 10 ). On the fine as well as on the coarse grained scale, however, quantum and classical transmission probabilities 
differ dramatically. 
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